Interaction-controlled transport of an ultracold Fermi gas 



O 

o 

(N 
o 

Q 
in 



Niels Strohmaier^, Yosuke Takasu^'^, Kenneth Giinter^, Robert Jordens^, 

Michael Kohl^''^, Henning Moritz^'*, and Tilman Esslinger^ 
^Institute for Quantum Electronics, ETH Zurich, 8093 Zurich, Switzerland 
^Department for Electronic Science and Engineering, Kyoto University, Kyoto 615-8510, Japan 



Cavendish Laboratory, 



University of Cambridge, Cambridge CBS OHE, United Kingdom 
(Dated: February 1, 2008) 



We explore the transport properties of an interacting Fermi gas in a three-dimensional optical 
lattice. The center of mass dynamics of the atoms after a sudden displacement of the trap minimum is 
monitored for different interaction strengths and lattice fillings. With increasingly strong attractive 
interactions the weakly damped oscillation, observed for the non-interacting case, turns into a slow 
relaxational drift. Tuning the interaction strength during the evolution allows us to dynamically 
control the transport behavior. Strong attraction between the atoms leads to the formation of local 
pairs with a reduced tunneling rate. The interpretation in terms of pair formation is supported 
by a measurement of the number of doubly occupied lattice sites. This quantity also allows us to 
determine the temperature of the non-interacting gas in the lattice to be as low as (27 ± 2)% of the 
Fermi temperature. 

PACS numbers: 05.60.Gg, 03.75.Ss, 71.10.Fd 
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The study of conductivity in solids has led to the dis- 
covery of fundamental phenomena in condensed mat- 
ter physics and to a wealth of knowledge on electronic 
phases. Intriguing quantum many-body phenomena such 
as superconductivity and the quantum Hall effect mani- 
fest themselves in their characteristic electronic trans- 
port properties. Moreover, the ability to manipulate con- 
ductivity has found numerous applications in technology, 
most prominently in semiconductors, where transport is 
controlled by electric biasing. Other ways to modify the 
conductivity in a material include adjusting temperature, 
pressure or magnetic field. 

A gas of ultracold fermionic atoms exposed to the po- 
tential of an optical lattice offers a new approach to study 
and control transport while providing a direct link to 
fundamental models in condensed matter physics. A pe- 
riodic potential of simple cubic symmetry is generated 
by three mutually perpendicular laser standing waves re- 
producing the potential experienced by electrons in the 
crystal structure of a solid. Prepared in two different 
spin states, fermionic atoms mimic spin-up and spin- 
down electrons. A unique feature of the atomic system 
is that the strength of the coUisional interaction between 
the two components can be directly tuned using a Fesh- 
bach resonance [J, [3] ■ While this property has been used 
to study fermionic superfluidity in the strongly interact- 
ing regime (e.g. d, 01), it has so far not been apphed 
to investigate transport phenomena in optical lattices. 
In previous experiments the transport of non-interacting 
fermionic atoms and the effect of a bosonic admixture 
mediating interactions were studied in one-dimensional 
optical lattices [EHQ- Furthermore, the dynamics of 
Bose gases in a three-dimensional optical lattice was in- 
vestigated experimentally and theoretically [^, 0, • 

In this letter we study the transport properties of 
a two-component ^"K cloud trapped inside a three- 
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FIG. 1: (Color online) Dynamic control of transport by tuning 
the coUisional interaction. The graph shows the center of mass 
motion of a two-component Fermi gas of (2.9±0.3) x 10* atoms 
in a lattice of 5 Er depth. At f = 0, the equilibrium position 
of the underlying harmonic trap is displaced vertically. After 
9 ms of evolution without interaction (•), the magnetic field 
is changed linearly in 1 ms so that the interaction is strongly 
attractive (o) for 10 ms. Then the magnetic field is changed 
to its original value again within 0.5 ms (•). The two non- 
interacting cases • are each fit by a damped cosine and an 
offset, while the evolution with attractive interaction o is fit 
by an exponential decay and an offset. Error bars denote the 
standard deviation of at least 4 measurements. 



dimensional optical lattice with underlying harmonic 
confinement. We monitor the center of mass motion 
of the atomic cloud after a sudden displacement of the 
trap minimum. The regimes of vanishing, weakly attrac- 
tive and strongly attractive interactions are accessed by 
exploiting a Feshbach resonance to tune the scattering 
length for low energy collisions between the two atomic 
components. The atom number is adjusted so that at 
the trap center the lowest energy band is either filled 
or half-filled. For these parameters the system can be 
regarded as a realization of the attractive single-band 
Fermi-Hubbard model llj with additional harmonic con- 
finement. A Mott insulating phase of pairs as discussed 
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in the context of the multi-band Hubbard model 0, [l3| 
is not expected. 

The strong influence of the interactions on the trans- 
port is illustrated in Fig.[TJ An atomic cloud is prepared 
in the optical lattice at half filling and brought into non- 
equilibrium by displacing the trap minimum. The ini- 
tially non-interacting cloud performs a weakly damped 
oscillatory motion in the confining potential. By tem- 
porarily switching on the attractive interaction, a con- 
trolled interruption of this oscillation is achieved. 

Our experimental setup that is used to produce quan- 
tum degenerate Fermi gases is described in detail in pre- 
vious work 13|. In brief, we prepare a cloud of 



atoms in an equal mixture of the hyperfine substates 
|F = 9/2, mF = -9/2) and \F = 9/2, mp = -7/2) in 
a crossed-beam optical dipole trap operating at a wave- 
length of 826 nm. After evaporative cooling we obtain 
4 X 10'* (3 X 10^) atoms at temperatures below T/Tp — 
0.20 (0.25) in the dipole trap with final trapping frequen- 
cies of {ujx,ujy,uJz) = 27r X (35, 23, 120) Hz, where Tp is 
the Fermi temperature. Next, the degenerate Fermi gas is 
subjected to the additional periodic potential of a three- 
dimensional optical lattice with a depth of 5 En. The 
recoil energy is given by Er ~ /i^/(2mA^), with h being 
the Planck constant, m the atomic mass and A = 1064 nm 
the wavelength of the lattice beams. The lattice is formed 
by three retro-reflected laser beams with circular pro- 
files having 1/e^ radii along the (x, y, z)-directions of 
(160, 180, 160) fim at the positions of the atoms and a 
mutual frequency difference of several 10 MHz. To load 
the atoms into the lowest Bloch band of the optical lat- 
tice we increase the intensity of the lattice beams using 
a spline ramp with a duration of 100 ms at a scattering 
length of a = 50 oq, where oq is the Bohr radius. 

The gas is brought into a non-equilibrium position by 
increasing the beam intensities of the underlying dipole 
trap, which shifts the trap minimum by up to 2.5 fim 
in the vertical z-direction. Since this displacement is 
smaller than our imaging resolution, we map the center 
of mass position of the atomic cloud to momentum space. 
For this purpose we switch off the optical lattice and let 
the cloud oscillate in the remaining harmonic dipole trap 
for a quarter period [§]. After free expansion, we obtain 
the momentum distribution of the cloud from absorption 
imaging, determine the center of mass momentum using 
a Gaussian flt and infer the original displacement z of 
the cloud in the trap. Oscillations of the cloud size are 
not observed since the horizontal and vertical trapping 
frequencies are only increased by about 4% and 12%, re- 
spectively. Also, Bloch oscillations can be neglected for 
our small displacement since even for large fillings only 
few atoms gain sufficient energy to reach the band edge 
The energy deposited in the system by the trap 
displacement is estimated to increase the temperature in 
the lattice by an amount of 0.05 Tp. 

Variation of the magnetic bias field in the vicinity of 



the Feshbach resonance at 202.1 G Q allows us to tune 
the coUisional interaction between the two components of 
the Fermi gas. Prior to the displacement of the trap the 
magnetic field is gradually ramped to final values between 
210 G and 202.95 G within 50 ms, yielding an s-wave scat- 
tering length ranging from to — 1500ao. Usin g th e de- 
scription of a Hubbard model for cold atoms (isl Il6| , this 
corresponds to an effective interaction strength U/J be- 
tween and —24. Here U denotes the on-site interaction 
energy of two atoms in a different spin state, and J is the 
matrix element for nearest-neighbor tunneling, which has 
a value of J « ft- x 290 Hz for our lattice depth. 
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FIG. 2: (Color online) Evolution of the center of mass position 
for different interaction strengths and fillings. The circles o 
(•) denote samples in the half filled (band insulating) regime. 
For each data point the position of the cloud with and with- 
out displacement was compared to eliminate long term drifts. 
The error bars denote the standard deviation of at least 4 
measurements. 

The results of transport measurements for three 
interaction strengths and two different fillings are 
shown in Fig.O We fit the function z(t) = 
Zosc cos(27r/i) exp(— /5i)-|-2;exp exp(— Fi) + zo to the data. 
For the half filled case, (3.6 ± 0.4) x 10* atoms are pre- 
pared corresponding to 0.46 ± 0.05 atoms per lattice site 
and spin state at the center of a non-interacting cloud 
[l^. In the case of high filling, the samples contain 
(2.7 ±0.3) X 10^ atoms and form a large band insulating 
core. The filling factor of the Bloch band is position- 
dependent due to the inhomogeneous density profile of 
the trapped gas. In the following, we will discuss the 
cases of zero, moderate and strong interaction. 

In the non-interacting case with half filling we observe 
damped dipole oscillations (Fig.[2Ka) o). This damping 
of the center of mass motion can be attributed to the 
fact that the fermions in different quasi-momentum states 
possess different effective masses, resulting in a spectrum 
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of oscillation frequencies. Furthermore, the total trap- 
ping potential is slightly anharmonic, which causes a de- 
phasing also observed in the pure dipole trap. The sys- 
tem with high filling (Fig.[5fa) •) is characterized by a 
very slow relaxation towards the equilibrium position: 
The band insulating core suppresses center of mass mo- 
tion and a large number of atoms occupy localized states 
BEIli. These single particle eigenstates exist at a dis- 
tance z\oc from the center of the trap where the potential 
energy due to the harmonic confinement is larger than the 
bandwidth, i. e. mu1zy^^/2 > 4J. Consequently, the mo- 
tion through the center is energetically prohibited, how- 
ever, the atoms can still oscillate within the outer regions 
of the cloud. Even in the half filled case a small fraction 
of atoms is localized, which explains the small offset ob- 
served in the center of mass position after the decay of 
the oscillations (Fig.[2Ua) o). 

For moderate attractive interaction and half filling, the 
damping of the dipole oscillations becomes more pro- 
nounced (Fig.[2Ib) o). The damping rate (3 increases from 
(80 ± 17) Hz in the non-interacting case to (140 ± 37) 
Hz at U/J = —2.4. As the interaction strength is in- 
creased beyond U/J< —3.5, the oscillations vanish en- 
tirely. The sample with high filling (Fig.[2l^b) •) relaxes 
faster towards equilibrium than in the non-interacting 
case, which can be attributed to umklapp processes [20| . 

In the strongly interacting case, a very slow relaxation 
is observed for both fillings [Fig. [21(c)]. The transport in 
this regime is governed by the dynamics of local fermionic 
pairs. In the limit of low atomic densities bound pairs 
form for U/J < -7.9 [llllll,!!!]. These pairs tunnel to 
adjacent sites via a second order process with an ampli- 
tude Jeff = 2J^/[/. This effective tunneling is obtained 
by considering a ground state where all atoms form pairs 
and by treating the tunneling term proportional to J as a 
perturbation in the Hubbard Hamiltonian ll|. Accord- 
ingly, the tunneling rate of pairs is reduced with increas- 
ing interaction as compared to bare atoms. Besides, the 
energy offset between neighboring sites due to the har- 
monic confinement reduces the tunneling probability. For 
these reasons we expect the relaxation time to become 
longer for stronger interactions. This is supported by the 
data in Fig.[3l which shows a clear decrease of the relax- 
ation rate F for growing attractive interaction. The data 
is well fit by the empirical power law T/J oc {U/ J)~^'^^ . 
A quantitative understanding of this behaviour is chal- 
lenging due to the coexistence of bare atoms and local 
pairs which act as hardcore bosons in the lattice. 

Further insight into the physics of local pairs is gained 
by probing the double occupancy in the lattice for var- 
ious interaction strengths without displacing the trap. 
For this purpose, we prepare the system at half filling, 
as before, and set the desired interaction within 50 ms by 
changing the value of the magnetic field. Then the lattice 
depth is abruptly increased from 5 En to 30 En in order 
to suppress further tunnelling. By subsequently ramping 
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FIG. 3: Relaxation rate F as a function of interaction 
strength. The data points and error bars are fit results to 
center of mass evolutions of (3.8 ± 0.4) x 10* atoms. The 
empirical power law hT = CJ{J/UY with u = 1.61 ± 0.03 
and C = 2.95 ± 0.17 fits the data well. The inset shows the 
dependence of the relaxation rate on the atom number. 



the magnetic field from 203.26 G to 201.23 G within 5 ms, 
weakly bound Feshbach molecules are formed on those 
sites which are doubly occupied |2l[. We determine the 
number of atoms remaining after the molecule formation 
and compare it with the atom number which is obtained 
after dissociation of the molecules by applying the inverse 
magnetic field ramp. This yields the molecular fraction 
displayed in Fig.|4l showing a strong dependence on the 
interaction strength: While for the non-interacting sys- 
tem the detected fraction is 18%, it increases up to 60% 
for strongly attractive interactions. 

For the non-interacting gas the double occupancy in 
the lattice is solely determined by the number of trapped 



atoms and their temperature [19|, l21|, |23| . The detected 
fraction of 18% is consistent with the temperature in the 
lattice of (0.27±0.02) Tp, which we determined in a sepa- 
rate measurement with 2.7 x 10^ atoms yielding a molec- 
ular fraction of (45.7 ± 2.4)%. This temperature, even 
though measured in an ideal gas, suggests that the gas 
remains above the critical temperature for superfluidity 
111 , |24| . |25| also for strong interactions. Numerical cal- 
culations for an homogeneous interacting system show a 
considerable temperature dependence of the double oc- 
cupancy [2J|. We therefore expect that the temperature 
of the interacting gas can be deduced from the measured 
double occupancy. 

The increase of the molecular fraction with rising at- 
tractive interactions provides strong evidence for the for- 
mation of local pairs. In accordance with numerical cal- 
culations for the attractive Hubbard model at finite tem- 
perature 24, [2^ the number of doubly occupied sites 
increases already for weak attraction, i. e. in a regime 
where no bound state exists in the two-body problem 
{U/J > —7.9). Pair formation in the many-body regime 
is expected to start at a value of J7/J « —2 [2^. 

For strong attractive interactions, i.e. U/J < —7.9, 
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FIG. 4: The fraction of molecules formed in the optical lattice 
increases with attractive interaction, demonstrating a higher 
number of doubly occupied lattice sites. Error bars denote 
the statistical errors of at least 4 measurements. The line 
serves as a guide to the eye. 

the number of doubly occupied sites saturates. This is 
in accordance with the fact that the pairs are well local- 
ized on single lattice sites and can be regarded as hard- 
core bosons. An increase in molecular fraction due to an 
attraction-induced shrinking of the cloud, which would 
result in a higher average density, is not substantiated 
by the following measurements: When tuning the inter- 
action strengths we could not detect a change in the size 
of the trapped atom cloud with our measurement accu- 
racy of 10%. Furthermore, the same increase in molecular 
fraction is found if the attractive interaction is turned on 
within only one tunneling time. This demonstrates that 
we observe local pairing rather than a redistribution of 
the trapped atoms on a larger scale. 

In conclusion, we have found that the transport of an 
attractively interacting Fermi gas in a 3D optical lattice 
is strongly influenced by the formation of local pairs. In 
the future, studying the oscillation frequency below the 
superfluid transition temperature could serve to charac- 
terize the BCS-BEC crossover [2^. Extending our stud- 
ies to the repulsive Fermi-Hubbard model may provide 
a tool to identify quantum phases such as the fermionic 
Mott insulator [22] . 
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